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UNIQUENESS OF NONNEGATIVE WEAK SOLUTION TO u p < (-A f-u ON R w 


YUZHAO WANG AND JIE XIAO 

Abstract. This note shows that under (p, a, N) e (1, oo)x(0,2)xZ + the fractional order differential 
inequality 

(t) U p < (-A)?M in R* 

has the property that if N < a then a nonnegative solution to (f) is unique, and if /V > a then 
the uniqueness of a nonnegative weak solution to (f) occurs when and only when p < N/(N - a), 
thereby innovatively generalizing Gidas-Spruck’s result for u p + Am < 0 in M. N discovered in UTll . 
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1. Statement of Theorem 

In their 1981 paper [Hill . B. Gidas and J. Spruck discovers that the Laplacian differential in¬ 
equality on Euclidean A-space R w with A e Z+ = {1,2,3,...}: 

(t) u p <-Au under pe(l,oo) 

obeys the following law: if N < 2 then a nonnegative weak solution to (£) is unique; and if 
N > 2 then the uniqueness of a nonnegative weak solution to ($) happens when and only when 
p < N/(N - 2). 

The fact that this result is extendable to more general differential inequalities in R' v and even 
Riemannian manifolds has been observed now by several authors with a variety of arguments; see 

u a m tni m m m \n\ m m m m m . 

In this note, we will establish an analogue of Gidas-Spruck’s discovery for fractional order 
Laplacian inequality. To state our main result, let us recall some notations and facts on fractional 
Laplacian. 
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For 0 < a < 2 let (-A)" be the ^ (fractional order) Laplacian which is determined via the 
Fourier transformation (•): if / satisfies 


( 1 ) 

then 


f 

Jr n 


\m\ 


(1 + |x|) w+ “ 


dx < oo 


(2) (-A)?/(£) = \fff(f) = l^r f e~^f(x)dx V 

J R" 

alternatively, one has 

(3) (~A)%f(x) = C N<a p.v. f ^ ~ V 

j R « |jc - y| 


with Cjv, a being a normalization constant. The fractional Laplace operators arise in stochastic 
theory as the operators associated with symmetric a-stable Levy processes. Applications also can 
be found in financial mathematics, such as the pricing model for American options, we refer the 
reader to [j2][8j] and reference therein. 


Theorem 1.1. For (p, a, N) e (1, oo) x (0,2) x Z+, the fractional differential inequality 
(4) u p < (-A)5« in R w , 

enjoys the following property: 

(i) If N < a, then a nonnegative weak solution to © is unique; 

(ii) IfN > a, then the uniqueness of a nonnegative weak solution of © occurs when and only 
when p < N/(N - a). 

Remark 1. The notion of a nonnegative weak solution to © is deferred to Definition l2.21 

Notice that the case a = 2 of Theorem 11.11 agrees with Gidas-Spruck’s uniqueness result stated 
above. To prove Theorem [Ll] we utilize Section 2 to introduce three new results of independent 
interest: the first is Lemma 12.11 which, as a local version of ©, is established by Caffarelli- 
Silvestre’s cr-extension in ©; the second is Lemma [231 that lifts L P (W) to a weighted function 
space in the upper-half-space via the o'-cxtcnsion; and the third is Lemma [2~4l which exists as a 
mixed trace inequality (inspired more or less by If22l0 for the o'-cxtcnsion of a nonnegative weak 
solution to ©. With the help of those three lemmas, Section 3 is designed to validate Theorem 
11.11 - in particular - Theorem 1 1 . 1 l i ) and the “when” part of Theorem 1 1 . 1 l ii ) are demonstrated by 
showing “w > 0 => u — 0”, and the “only when” part of Theorem I Lill i) is checked through 
constructing a positive solution based on the fundamental solution to (-A )%u = 0 in R N . 

Notation. In what follows, U < V means that there is a constant C > 0 such that U < CV. 
Moreover, U * V stands for U < V < U. 


2. Lemmas and Their Proofs 

2.1. A reduction for u p < (-Apu. Such a reduction comes from another look at Caffarelli- 
Silvestre’s extension technique (cf. ©) that reduces the fractional Laplacian to a local problem 
through bringing one more variable into play. To be more precise, the fractional Laplacian can be 
characterized as a Dirichlet-Neumann operator for an appropriate differential equation of diver¬ 
gence form: if a e (0,2), u = u(x) is a function defined in R N , and to = to(x, y) is a solution to the 
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boundary value problem 

(-diviy^Vco) = 0 in R^ +1 = R N x (0, oo); 

0) = u(x) V x eR n , 

then there is a constant C a depending on a such that 

(6) lim y 1- "— (x,y) = —C a (—A)%u(x) V x e R w . 

y-> o + oy 

Observe the divergence form in © is not only the Euler-Lagrange equation for the functional 

J(u) = ff |Vcu| 2 y 1- " dxdy 
JJ R^ +1 

with V being the gradient in (x, y). but also it can be rewritten as 

(7) A x (j) + (1 - a)y~' OJy + (Oyy = 0, 


which may be viewed as a harmonic function in N + 2 - a dimensions. This observation leads us 
to the fundamental solution to ©: 


7r w+2-ff)/2r((jy - a)/4) J A7 

r(x, y) =-—- under N -a > l. 


N-a 

,123 — 


(W 2 + Ivl 2 )' 

where T(-) is the standard gamma function. Furthermore, upon setting 

, -!« OT v ^ W(N + a)/2)\ y 

P(x, y) = -y d y T(x, y) = - 

l ^ /2 n«/2) / (\ x \2 + I 


(|.v| 2 + |y| 2 )" 


be the Poisson kernel, we find that 

(8) co(x, y) = f 

J r n 


P(x - z, y)u(z ) dz. 


solves ©. Due to ®, the solution a>(x,y) of © is also called the cr-extension of u(x). 

One major advantage with the use of representing the fractional Laplacian as a divergence- 
form differential equation ©© instead of a nonlocal operator © or © is that one can local i / e 
the related problems. This method has a great deal of applications in equations involving the 
fractional Laplacian; see © IS HI and their references. Our current concern is the inequality © 
(but not its equality) that is also nicely connected to the extension method in [©. After extending 
the inequality to JEl^ +1 , it suffices to deal with a local equation © by using an inequality as the 
boundary condition. In view of © and ©, the inequality © can be rewritten as 


(9) 


div(y l ~ a Wco) = 0 V (x,y) e 


sAT+l. 


dco 


lim y l ~ a —(x,y) + u p (x) < 0 V x e 
y-> 0 + oy 


»v 


Naturally, co can be extended to 3c' v 11 via putting 


cb(x,y) 


I oj(x,y), 

| oj(x, -y), 


V x e R n & y > 0; 

V x g & y < 0. 
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Lemma 2.1. Given (p, a) E (1, oo) x (0, 2). Suppose that u is a nonnegative solution to ©, co is 
its a-extension, and V stands for the full gradient operator in (x,y) e R w+1 . If \y\ x ~ a \V co{x,y)\ G 
L 2 (R W+1 ), then for any nonnegative continuous function h satisfying |y| l_ °'|V/i(v, y)| G L 2 (R W+1 ), 
one has 


( 10 ) 


f 

Jr n 


u p (x)h(x, 0) dx 


*rr 

JJm.n+1 


Vcb(x,y)-SIh(x,y)\y^ 01 dxdy. 


Proof Without losing of generality, we may assume that h supports in the origin-centered Eu¬ 
clidean ball B r := B(0,R ) in R^ -1-1 with radius R > 0. For any s > 0, we utilize div(\y\ ] ~""V co) = 0 
and the Stokes formula to obtain 


ff. 


Vco(x,y) ■ Vh(x,y)\y\ a dxdy 

= // ♦// 

JJb«\(M< £ } JJb , 

-ff 

JJb, 

-ff 


Va>(x,y) ■ V/i(jc,y)|y| 1 a dxdy 


OJB R n{[v|<£) 

div(\y \ 1 ~ a hVco)dxdy + 
deb 

s x ~ a h(x, s)—(x, s)dx + 
B«n(|y|= e | oy 


ff 


Vco(x,y) ■ V/?(A,y)|y| 1 a dxdy 


ff 

JJb, 


Vco(x,y) ■ V/i(x,y)|y| 1 01 dxdy. 




Upon noticing that the last integral goes to 0 as s —> 0 + because of lyl 1 “|Va>(;c,y)| G Lf oc (R N+1 ) 
and that the inequality in © gives 

u p (x) < - lim s 1 <r -f—(x, s), 

£ —>o + oy 

we get (fTOl) . □ 

Remark 2. As a matter of fact, when h(x,y) is the a-extension of h(x, 0), the inequality (fTOl) is 
equivalent to 


f 

Jr n 


u p (x)h(x, 0 )dx < f \f\ a u{f)h(f, 0 )d£, 


which implies u e H a/2 P\L p+1 immediately. One way to verify this is to show that the correspond¬ 
ing energy functionals coincide through the energy functional formula in 0: 


(ID 


ff 

Jok'H 


\S7co(x,y)\ 2 y x “ dxdy 


f 

Jr n 


mmvdf. 


where co is the a-extension of a. A similar argument gives 


ff 

JJrH 


Vco(x,y)-Vw(x,y)y x a dxdy 


-f 

J r n 




where v g H al2 and w is its a-extension. 


More importantly, Lemma 12.11 suggests us to adopt the local version of a nonnegative weak 
solution to © via the a-extension. 

Definition 2.2. We say that u on R N is a nonnegative weak solution of © provided that u is a 
nonnegative function and its extension co satisfies both |y| l_ ' K [Vto(x,y)l e L 2 (R n+1 ) and 


(12) 


f u p h\ y=0 dx < ff Vco ■ Vh lyl 1 01 dxdy, 

J RN JJ Ri' +1 
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for any compactly supported nonnegative function h satisfying |y|' a \Vh(x,y)\ e L 2 (R N+1 ). 


2.2. An extension for L p { R"). This is presented in the following result. 


Lemma 2.3. For ( p , a) e (1, oo) x (0, 2) and u e L P (R N ), let to(x,y) be given by ®. Then 

(13) ^JJ ^ Ia>(x,y)l " y la dxdyj < IMIzar")- 

Proof. The argument is split into two parts. 

Part 1: 2N/(N + 2 - a) < p < oo. This follows from the well-known complex interpolation 
theorem. So it suffices to prove the endpoint cases, which are (oo, oo)-estimate and (2, N l ^_ ir )- 
estimate. 

Let us first consider the (oo, oo)-estimate. According to ®, one has 


oj(x, y ) = 


T((/V + a)/2) \ r y°u(f) 
7r N P-r(a/2) ) Jrn (fx - fj 2 + y 2 )' 


d£- 


It is easy to get 


f 

Jr n 


y 


(\x\ 2 +y 2 )- 
and then by Young’s inequality 


dx 


= (■—*, 

Jr w (1 + |x| 2 )" 


dx < oo, 


ll^(-Ty)llL”(R'v +i ) ^ II^WIIl“(r' v )- 


This in turn implies 


(14) \\to(x,y)\\ L7y 

(R ¥ +l ,y 1 - a dxdy) < \\u(x)\\ L ™. 

Next, we make the //-estimate with p = ■ After taking the Fourier transform, d7]) becomes 

-\£\ 2 cb(f,y) + (1 - a)~ ] ytOy(£,y) + tb yy (£,y) = 0 , 

which is an ordinary differential equation for each £. Suppose now that 0 : [0, oo) —» R is the 
minimizer of the functional 

(|0| 2 + |0'| 2 )y 1_ff dy subject to 0(0) = 1. 
o 

Then 0 solves the following equation (cf. Il4l). 

-0(y) + + (pyy(y) = 0; 

• 0(0) = l; 

lim v ^oo 0(y) = 0. 

Note that 


oMf, y) = u(£)f(\£\y). 
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Thus the energy functional of co is 


ff 

JJ R* 


\u(x,y)\ 2 y l a dxdy = f f 
Jr. n Jo 


m,y)\y- a dyd£ 


\mmm\y- a dyd{ 


(15) 


I^OOlV a dyd% 


- f r 

= r i^r- 2 ia(^)i 2 r 

Jrw Jo 

<J((f))\\u\\ 2 a t 

H 2 

<J(cf))\\u\\ 2 2 ,v , 

LN+'l-a 

where we have used the fractional Sobolev imbedding for | - 1 < 0 in the last step. 

The desired inequality follows from putting the above two estimates together and using IfTOl 
Theorem 1.3.4] with 

(X,dp) = (R N+] ,y ] ~ a dxdy); (Kdm) = (R N ,dx); (q 0 , Po,Qi, Pi) = ( 2 , -, 

\ N + 2 - a 

and T : u h» co being the extension operator. 

Part 2: 1 < p < 2N/(N + 2 - a). This follows from the well-known real interpolation theorem. 
So it suffices to check the weak-type estimate at the endpoint p = 1. 

Given a space X and a measurable function / : X —> C, we define the distribution function 
Af : R+ —» [0, + 00 ] by the formula 

A f (t)=p({xeX:\f(x)\>t}). 

Recall the weak L p -norm W/Wlp^cx) by 

:= 11 ^/( 011 ^ 0 ^)- 

We say that T is of strong-type ( p , q) if 

\\Tf\\mx,dr) ~ \\f\\LP(Y,dm), 

and of weak-type (p, q) if the above inequality holds whenever L q is replaced by ZJ °°. 

Suppose now dp = y l ~ a dxdy. Then the weak //“-norm has an equivalent counterpart 


(16) 


sup p(Ey 


£ 


f dp 


where the supremum is taken over //-measurable set E c Wp . If r 0 = ^ - , then for Tu = co we 

have 


\\Tu\\ 


L r 0' c 


: SUp p(E) r o 
E 


£ 


Tudp 


For the right side of the last equivalence, we use Fubini’s theorem and (fl6l) to achieve 


p{E) '0 


£ 


Tudp 


<p{E) 


■*££ 

Je Jr n 


P(x - z,y)\u(z)\dz dp 


= |sup//(£ - z) 
<\\P\\ L r 0 - \\u\y. 


y ■ 


P(x, y) dp 
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where z = (z, 0) g R+ +1 . 

Now for the weak (r 0 , l)-estimate of T, it is sufficient to bound 


namely, 

(17) 


\\P\\u»-°°{X4n) for P(x,y) 


m(N + a)/2) \ y- N 
{ n N/2 T(a/2) j (|jc/y| 2 + 1)^ ’ 


H p H^(^) = supd>({(x,y) g Rf 1 : P(x,y) > A}). 


Noticing 


P(x,y) 


m(N + a)/2)\ v" 

1 F(a/2) 1 ( 1 ^+/)*? 


we find that the right side of ([T71) can be controlled by 


where X = (x,y). 


mx,y) 6 Rf 1 : 


P{x,y) > A}) <p ({X G Rf- 1 : |X|- W > d}) 


< 


f , 

dp 

J\X\<A N 


(f ■ 

dx)( f dy ) 

W|x|:ST N 

/ \ Jy<A * ) 

A ro , 



thereby getting 


ll^llz>“(X,d/i) ^ 1, 


which implies that T is of the weak (r 0 , l)-estimate. Therefore the strong-type estimate desired 
in Part 2 follows from the strong (2, v j^ v ^ -estimate proved in Part 1 and Marcienkiewicz’s real 
interpolation method IfTOl Theorem 1.3.2], □ 


Remark 3. In fact, if we interpolate between the weak (r 0 , l)-estimate and the strong (oo,oo)- 
estimate, the strong-type estimate for 1 < p < oo can be obtained. However, the reason that we 
split the argument into two parts is that not only the proof of the first part is more direct, but also 
the exploited energy structure is a very powerful tool and has its own interest. 

2.3. A mixed trace estimation for u p < (-A )?u. Such a type of estimate is given below. 


Lemma 2.4. Let ( p,a ) G (l,oo) x (0,2) and ip be a smooth function in R+ +1 with the compact 
support supp(if) such that 0 < tp < 1 and p = 1 in an nonempty open subset of supp(<p). If u is a 
nonnegative weak solution to © and to is its a-extension, then for any 0 < t 1 and s » 1 there 
is a constant C v ,v > 0 depending on s and N such that 
( 18 ) 


/ 

Jr n 


*P lv=0 


r dx < C 




|V^| 


2p'-2t' 

p'-l+t-t 1 


y l ~ a dxdy 


where 


/ 

P 


(N + 2 - a)p 


& 


(N+ 2- a)t 


N 


N 
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Proof. Since u is nonnegative, so is co. For constructing a proper test function, we introduce a 
small number 

0 < 8 <sc 1, 

and let 

(19) u) S = (j) + 6 & i//(x,y) = 

where 0 < t <c 1, s » l.A simple calculation derives 

(20) Vt/r = -fo>J 1- y Va> + sa)f(p s ^ l Vip. 

Using the compactness of supply), we obtain f e L P+1 (R N+1 ). And since co s is uniformly away 
from 0, a>g l is uniformly bounded from above, whence getting that for the fixed t > 0 and .v > 1 

\Vf{x,y)\ 2 y l ~ a dxdy e L 2 (< +1 ). 

The foregoing analysis tells us that such a f can be chosen a legal test function h for (fT2l) . So, 
by (fT2l) and (12(71) we obtain 

( 21 ) 

t ff aj-'- t ip s \Vajfy ] -°dxdy+ f <p s \ y=0 (u + d) p " ? dx < s ff tof ip s ~\V oj • Vip)y l ~ a dxdy. 

JJr " +1 Jr w ' JjRy +1 

By the Holder inequality and then Young inequality, the right hand of (I2TI) can be estimated as 
follows: 


-JL 


Va> ■ V^y 1 a dxdy 


(tu) s l '(p s )HV(o-V(p)st L/Z u>' 2 ip 2 1 y 1 a dxdy 


1/2 If . l-l \ - a . 


< U 


ff a)- s l - , (p s '\Vco\ 2 y l - a dxdy) (rf 1 ff aYf 1 ip s ~ 2 \V ip\ 2 y l ~ a dxdy\ 
JJr ^ +1 / \ JJrJ ,+1 / 


JL 


ujg 1 VlVwlV “dxdy + 


^ ff 

2 JJr" 


ojg y 2 |Vy?| 2 y l a dxdy 


Bringing this last inequality into (12TI) we obtain 

( 22 ) 

*- ff (o- s 1 - t ip s \Vu)\ 2 y l -° , dxdy+ f <p s \ y=0 (u + 6)^ dx < ff co]-'w s ~ 2 \Vcp\ 2 y 1 - a dxdy. 

2 JJrJ'+i JrN 2 JJ R «+1 

By the Holder inequality, the right side of (l22l) can be estimated as follows: 

2 - yr‘ ff aj 1 s -y- 2 IV<fl 2 y 1 - a dxdy 

JJ R ^ +1 

( 23 ) = ff oj f y S (2 ~ 1 .s' 2 r 1 ,s | V(/j| 2 ) y 1 ~ a dxdy 

JJ< +1 v ' 

where 


P i 
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Note that for p > 1, we have 


P' > 


N + 2-a 
N 


(, s 2 t 1 


And for s » 1 and 0 < t <sc 1, we have 

s— 2 — 7-77 ^ / 1 

p~ f < 1 . 

So, combining (l22l) and (1231) we finally obtain 


f <p s \y=o(u + 6) p 1 dx 

J r n 

JIT a/- i '<p s y 1 - a dxdy) P I ff t~^ff \V(p\7^ y'~ a dxdy 
\JJ r^ +1 / \JJs^ +1 JJpy +i 


p'-l+t-t' 

7=? 


Letting 6 —> 0, we obtain (fl8l) . thereby completing the proof. □ 

3. Verification of Theorem 

Proof of d^e case N < a or N > a & p < N/(N - a). Assume that u is a nonnegative weak solu¬ 
tion to ©. Then u satisfies (fl2l) . If to is the ar-extension of u, then Lemma \2A\ holds. 

In what follows, we estimate the second factor of right side of (fl8l) by selecting a series of 
appropriate test functions. To be more explicit, for a large number R > 0 let us consider the 
function 


(24) 

where 


<p(X) = 


L 


\X\ < R, 


|f| , \X\ >R, 


X = (x,y) e Rf 1 & |A| = yj\x\ 2 +y 2 . 

Since R is big enough, t = (InT?) -1 is sufficiently small. For any n e Z + define a cutoff function rj n 
by 

T, 0 < |X| < nR; 


(25) 

and then consider the function 


r 1n (X) = {2-§, nR<\X\<2nR- 


0, 


\X\ > 2nR, 


<P«Q 0 = <P(X)rj n (X) 

so that ip n (X) approaches to tp(X) from below as n —> 00 . Notice that 

V «fn = </? V l] n + TJ n V<fi. 

Thus, for any a > 2 one has 

(26) 1 vip n \ a < i^v/7„r + \n,y^\ a . 

Now we are ready to estimate the left side of (fT8l) by choosing <p n to be the test functions. More 
precisely, we need to deal with the integral 

| Xip u \ a y l ~ a dxdy. 


IrfjX) 


: ~~ff 
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From (1261) we have 
I n (a) < 


(27) 


ff r ] a n \V ( p\ a y l - a dxdy+ ff y‘\Vr ]n 
JJ Ri r+1 JJr" +1 


a 1X7^, \a^,l—a 


< 


ff 

JjR'i 


| V(p\ a y l ~ a dxdy + 


1 \b r 


f£ 


y dxdy 


<p a \^Tln\ a y X ~ a dxdy. 


( B2nR\B„ R )nWl 


For the second integral in (1271) . we use |V? 7 „| < (nR) 1 to gain 


ff 

JJ R" 


<p‘Wy l " dxdy <(nR)~ 


(28) 


<(nRY 


f 

J(L 


ip a \X\ l ~ a dX 


(B2hr\b „R)riS,+ 


(nR)'~ a \B(0,2nR)\ 


sup (f 

y(B 2 „R\B llR )nRy l ) 

d \ -at 

—j 0 nR) l ~ a (2nR) N+l 

N+2-a-a—at ryN+2-a-a 
~/7 A 

For the first integral in (1271) . we employ |Vyp| < RH\X\~ x '' to obtain that 


IL 


\V<p\ a y l ~ a dxdy < 


1 \b r 


(29) 


/ 

f 

Jr, 


\Xip\ a \X\ l ~ a dX 


r n+1 \b r 

Ra‘ t a r a-at r N + \-» dr 


R/2 

u ^ a j^~ a +N+2-a 


provided a + at>N + 2-a. Combining 


with (1271) yields 


I n (a) < n N + 2-a-a-at R N + 2-a-a + f a R -a + N + 2-a, 

and consequently, 

(30) I n (a) < n~ at + f under a > N + 2 - a. 

Here, it should be pointed out that I n {a) is uniformly bounded in R and n. 

Now, from (1181) and (l30l) we obtain a constant C > 0 depending only on 5 and a = such 
that 

f tpy-’ dx < C(n~ at + t a ) ( ff o/-'<p s n y l ~ a dxdyf ' . 

Jr n \JJ R*f +1 / 

Letting n —» oo, we further obtain 

(31) f tp s u p ~’ dx < t a ( ff a> p '~ l '<fi s y 1 ' a dxdy\ 

Jr n \JJrr +i I 

Since u e L p l (R N ), we conclude that m e L v (R^ +1 , dp) by Lemma [231 and that the 
integral of right side of OTT) is uniformly bounded in t. By letting t —> 0+, we discover 


/ 


(p s u p dx = lim f tp s u p 1 dx < lim C ( ff 
?_>0+ Jr^ f_ ^ 0+ \ J Jr" 

whence reaching w = 0 provided a + at>N + 2-a. 


a> p '~ t ' ip s y x ~ a dxdy ) = 0, 


1-r 

p'-t' 






UNIQUENESS OF NONNEGATIVE WEAK SOLUTION TO uf < (-A)?w ON R N 


11 


Note that 


a = 


2 <j>' - n 

p' -1 


So, choosing t to be sufficiently small one gets 


. . 2 p' , „„ 2 tf ^ 2 p' 

a + at = -1-2 1 -> 


P'~ 1 


p' - 1 p' - 1 


Thus we should guarantee 


2 P' 

p' - 1 


Recall that 


Thus we should have 


> N + 2 - a, i.e., (N - a)p' < N + 2 - a. 


, p(N + 2 - a) 
P = - N -' 


p(N + 2-a)(N-a) ^ AT , „ . p{N - a) ^ 

-< N + 2 - a, i.e., -—-< 1, 


N 


N 


as required. 


□ 


Proof of the case p > N/(N - a). Trivially, u = 0 is a solution to ©. But, a positive solution of 
© will be provided below. 

To do so, for 0 < s < N let I s be the Riesz potential determined by 

I s f = (-A )"*/. 


Actually, one has 


hfix) = 


T(N/2 - s/2) 


n N / 2 2 s T(s/2) 
Of course, the following fact is well-known: 


f 


f(y) 


& \x - y I 


N-s 


dy V x e 


»N 


(32) 


\x\ 


-N+s - 


= C N Jx\~ S & I s m = C N , S ^ nN l2 2sYis/2) j |x| ~ ./ (x) 


with c Ns being a constant. 

So, it is easy to find that the Riesz potential is the inverse of the fractional Laplacian in the sense 
of: 

4((-A)f/) = (-A )Hl a f) = f. 

By setting / = do (the Dirac point mass at the origin) in the last formula we get 


(-A)Hl a 6 0 ) = d 0 , 


whence seeing that 

, T s , \ (T(N/2-s/2)\ N 

(33) WAx) = i Mx) = y——— j W 

is the fundamental solution of (-Apu = 0. 

The main idea of the remaining argument is to perturb the fundamental solution (1331) properly 
to construct a positive solution to © when p < N/(N - a). However, the fundamental solution is 
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not a proper solution due to its singularity at the origin. So, in order to handle this issue, we need 
a smooth cutoff function p{x ) that is defined by 


r 1, \x\ < l; 

p(x ) = < smooth and radially decreasing, 1 < |x| < 2; 

[o, W > 2. 


It is easy to see that the function 
u(x) = p* W a (x) 


( UN/2 - 5/2) \ r p(y) 


XEl" 


is well-defined when 0 < a < N, and solves the a -fractional Laplacian equation, i.e., 

(-A)V*) = P(x) V x G R N . 


Thus, for 0 < 8 < N - a we perturb p * W ff (.r) at infinity to get 

us(x) := p * - -(x) = I" -——-dy V x e R N . 

I • \ n ~ a ~ 6 J R » lx - y\ N ~ a ~ 6 

We next claim that u& solves the inequality © under p > when |.r| is large. We observe 
u§(x) = I" -——-dy ~ ■——- V \x\ » 1, 

X\<i\x-y \ N - a - 5 

whence getting 

(34) u p Ax) ~ —7 -— V \x\ » 1. 

On the other hand, we use (l32l) to produce a constant c > 0 such that 

mo = c\d\- a ~ 6 p(a 

and thus 

(-A)T M<5 0 = c|f|-*p(£), 

which in turn implies 

2 1 

(-A ) 2 u s (x) = cp* 

Similarly, we have 

(35) (_A)i M|5 (*) * y | X | » 1. 

If p > jr~ r . then there exists a sufficiently small constant 6 such that (N - a - 8)p > N - 8. This 
inequality, along with (l34l) and (031) . implies 

ufj(x) < (-A)-u s (x) V \x\ » 1, 

thereby verifying the claim. 

Since we have showed that there is a positive R such that u s solves © with p > when 
|jc| > R, for the rest part |x| < R we need to exploit the scaling structure of the inequality. To be 
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more precise, observe that both u§ and (-A)'htg are continues and positive, and decay radially. So, 

letting 

W<« t 

- a -, 1 > , 

min(-A) T « <5 (x) 

M<« J 

we utilize the previous construction to deduce that ug solves the inequality 

u p < m{-A)~-u in R N . 

Finally, choosing u = m~u 6 we obtain a positive solution to © on ic' v , as required. □ 
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